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Abstract 

The  propagation  of  electromap;netic  pulses   around  the  earth  is 
investigated  analytically.      The  pulses  are  assumed  to  be  produced  by  a 
vertical  electric  or  map^netic  dipole.     The  earth  is  treated  as  a  homogeneous 
sphere  of  either  finite  or  infinite  conductivity  and  the  atmosphere  is 
assumed  to  be  homogeneous.      It  is  found  that  very  short  pulses  become 
lont^er  the  further  they  propagate,   in  addition  to  diminishing  in  amplitude. 
The  duration  of  a  pulse  which  is  initially  a  delta-function  increases  as  ^, 
where  &  is  the  angle  between  source  and  receiver.     The  results  are  represented 
as  products  of  several  factors,   which  we  call  the  amplitude  factor,   the  pulse-shape 
factor,   the  time-dependent  height-gain  factor  for  the  source  and  receiver,   and 
the  conductivity  factor.     Graphs  of  these  factors  and  of  the  pulse  shape  for 
several  cases  are   given. 
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1.  Introduction 

Suppose  an  electromagnetic  pulse  of  given  shape  and  amplitude  is 
emitted  by  a  source  on  or  near  the  ground.  We  wish  to  find  its  shape,  amplitude 
and  arrival  time  at  any  point  beyond  the  horizon  on  or  near  the  ground.  In 
particular  we  wish  to  determine  the  effects  of  the  ground  conductivity  along 
the  propagation  path,  Vte  will  assume  that  the  source  is  a  vertical  dipole  of 
either  electric  or  magnetic  type  and  that  the  earth  is  a  homogeneous  sphere. 
By  a  slight  modification  of  our  procedure,  we  could  also  treat  the  case  in 
which  the  ground  conductivity  varies  along  the  propagation  path. 

The  method  of  solution  consists  of  two  steps.  First  we  must  obtain  the 
Hertz  vector  which  represents  the  field  due  to  a  time  harmonic  or  periodic 
source.  Then  by  fburier  superposition  of  these  periodic  fields  we  can  obtain 
the  Hertz  vector  for  a  delta-function  source.  Other  pulse-type  sources  can  be 
treated  by  superposition  of  delta-function  fields. 

Since  the  field  due  to  a  periodic  source  has  been  determined  by  numerous 
authors,  our  first  step  is  simple — we  need  merely  copy  the  Hertz  vector  of  this 
field.  This  is  done  in  Section  2,  the  books  of  V.A,  Fock'-  ■'  and  H.  Bremmer'-  ^ 
being  used  as  sources.  In  Section  3  we  write  the  Fourier  integral  representing 
the  Hertz  vector  of  the  field  due  to  a  delta-function  source.  We  then  evaluate 
this  integral  by  the  saddle  point  method  and  obtain  a  simple  formula  as  the 
result  for  each  type  of  dipole.  In  Section  h   some  graphs  based  on  these  formulas 
and  a  discussion  of  the  results  are  given. 

We  wish  to  acknovjledge  our  indebtedness  to  Prof.  Bernard  Friedman  who 
previously  analyzed  the  field  due  to  a  pulsed  line  source  near  a  perfectly  con- 
ducting cylinder, 

2,  The  field  of  a  periodic  source 

Suppose  a  vertical  electric  dipole  is  located  at  the  point  r  ■»  p,  9  -  0 
of  a  polar  coordinate  system  with  the  origin  at  the  center  of  the  earth.  Then 
its  field  can  be  expressed  in  terms  of  a  Hertz  vector  having  only  a  radial 
corrponent  rU  (r,9,t),  by  the  fomulas 

(1)    E  =  V  X  V  X  (rUg) 

H  =  (e  ^  +  a-)  V  X  (?U^). 

Similarly,  the  field  of  a  vertical  magnetic  dipole  at  r  =  p,  9  =  0  can  be 
expressed  in  terms  of  a  radial  Hertz  vector  rU  (r,e,t)  by  the  formulas 


The  nvanbers    t       and   Z     are  complex,  lying  on  the  ray  arg   'T"  ti/3.     The    V^s 
of  smallest  absolute  value  are 

(16)  t^  -  0.808  e^"/3 

(17)  t°  -  1.856  e^^/3. 

For  a  highly  conducting  earth  cr,  is  large  so  6  is  large  and  6  is 
small.  Therefore  the  series  (lii)  is  useful  for  computing  T'  in  the  electric 
case  and  (1$)   is  useful  in  the  magnetic  case.  For  the  same  reason  (8)  is  more 
convenient  for  computing  f  in  the  electric  case,  and  (9)  is  in  the  magnetic  case, 

3.  The  field  of  a  delta-function  source 

The  Hertz  potential  U(r,©,t)  of  the  field  due  to  a  vertical  electric  or 
magnetic  dipole  having  any  time  dependence  can  be  obtained  by  Fourier  superposition 
of  the  periodic  potentials  u(r,9)e~    described  above.  Since  k  =  co/c,  where 
l/c  «  /efl  ,  we  can  integrate  with  respect  to  k  instead  of  «  and  obtain 


(18)    U(r,e,t)  -  ■§-     I   A(k)u(r,9)e-^''*  dk 


''-00 


In  (18)  the  Fourier  amplitude  A(k)   is  determined  by  the  time  dependence 
and  amplitude  of  the  source.     With  the  function  u(r,9)  defined  as  in  the  preceding 
section,   a  delta-function  source  has  the  amplitude  A(k)  -  1.      In  the  absence  cf  the 
earth,       the  Hertz  potential  of  the  field  due  to  such  a  dipole  would  be 

6(t  -  ?) 

o  '^       UtiR 


(19)         U     -  ° 


Here  R  =  (r  +  p  -  2pr  cos  B)    '     is  the  distance  from  the  source  to  the  receiver. 
Thus  with  A(k)  -  1,  (18)  will  yield  the  diffracted  pulse  due  to  the  incident 
pulsQ  (19). 

To  evaluate  (l8)  we  introduce  the  time  T  defined  by 

(20)    T  .  t  -  f  .  I  cos-^(|)  *f  cos-^(f). 

This  time  T  is  just  the  time  measured  from  the  arrival  of  the  diffracted  wave  front 
at  the  point  (r,9),   assuming  that  the  source  starts  at  t  =«  0,     We  now  insert  (8) 


$ 


or  (9)  into  (18)  and  set  A(k)  =  1. 

Let  us  first  consider  the  case  of  a  perfectly  conducting  earth.   In  this 
case  t   is  independent  of  k  and  5  is  either  infinity  or  zero.  Then  for  T  <  0 
the  contour  in  (18)  can  be  closed  in  the  upper  half-plane  to  yield  the  value 
zero  for  U.  This  is,  of  course,  to  be  expected  from  the  definition  of  T.  For 
T  >  0  we  evaluate  (18)  by  the  saddle  point  method.   In  the  electric  case 
we  use  (8)  in  (18),  set  6^  »  oo  and  obtain 
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In  the  magnetic  case  we  use  (9)  in  (18),  set  6=0  and  obtain 


(22)     U  (r,6,t)   -X 


co|r°l^/^ 


3      2^/2  2I/3  3I/2  „^l/2   (^^  ^)l/2 


(cT) 
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exp 
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(cT) 
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v(|r°l2^/3(^-l)}      v[lr°|2^/3(e(ga),^)j 
v'(-2^/3   |^0|^  V(-2^/3   |^0|) 


The  function  v  in  (21)  and  (22)  is  the  imaginary  part  of  the  /iry 
function,  v  =  ]jnw.  For  large  positive  z,  v(z)  is  asjTnptotic  to  ^  z~^'  exp[-  4  z   J 

while  for  large  negative  z,  v(z)  is  asymptotic  to  (-z)"  '  sin[»(-z)3''^  *  tI  *  ^" 
the  intermediate  region  the  values  of  v(z)  have  been  tabulated  by  V.  A.  FockW, 
Now  let  us  consider  the  case  of  a  finitely  conducting  earth.  In  the 


electric  case  we  will  use  for    t     the  first  two  temis  of  (lU)  with  6 


Upon 


inserting  (8)  into  (18),  with  this  value  of  "^  ,  and  putting  A(k)  ■=  1,  we  obtain 

an  expression  for  U  (r,6,t).  For  T  <  0  this  expression  is  zero,  as  before.  For 

T  >  0  we  again  use  the  saddle  point  method.  Since  6  is  a  function  of  k,  the 

determination  of  the  saddle  point  is  more  complicated  than  in  the  previous  cases. 

Therefore  we  neglect  the  term  -l/(2t°"  6  )  compared  with  t^   in  deteimining  the 

s   e  3 

saddle  point.   Then  we  obtain 


(23)     U^(r,&,t) 


'^P^t'^  na^/2(^.^  ^)l/2   •   21^00, 
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In  (23),   6"  is  the  value  of  6     at  the  stationaiy  point.     Its  absolute  value  is 


given  by 


(2li)      |5!|   » 
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Here  Tg  =  ©a/c  is  the  travel  time  for  the  pulse  to  reach  the  point  9  from  the 
source  at  &  =0. 

In  the  magnetic  case  we  use  the  first  two  tenns  in  (15)   for    f    with 

s 

6  =  6jjj.     Then  using  (9)  in  flS)  with  K^X)  =  1  we  obtain  an  expression  for 
U_(r,e,t).     As  before  this  expression  is  zero  for  T  <  0,  while  for  T  >  0  we 
obtain 


ce|r 


0|l/2 
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Here  6^  is  the  value  of  6  at  the  stationary  point  and  |6  |  is  given  by 


(26)    |6;|  - 


3CT  vl/2 


irri 


3'  /   ./e,      -^HcJ^licT)!:^ 

V2    i/^  ^V2 


U,  Discussion  of  results 

We  have  calculated  the  Hertz  vector  of  the  field  due  to  a  vertical  electric 
or  magnetic  dipole  with  a  delta  function  time  dependence.  For  a  perfectly  con- 
ducting earth  we  obtain  (21)  and  (22),  while  for  finite  conductivity  we  obtain 
(23)  and  (25).  Of  course  (23)  reduces  to  (21)  and  (25)  to  (22)  as  o-  -►oo. 
These  results  are  useful  at  points  beyond  the  horizon  and  therefore  only  the 
first  term  (s  =0)  in  each  result  need  be  retained.  This  is  so  because  the 
other  terms  have  relatively  rapidly  decaying  exponential  factors.  Consequently 
our  results  are  rather  simple  formulas. 

Let  us  examine  U^(r,9,t)  for  the  perfectly  conducting  case  when  both  the 
source  and  observation  point  are  on  the  ground  (p  =  a  and  r  =  a).  Then  (21) 
becomes 


(2?)  U^(a,9,t)  =  />g(e)  S(TAg^), 


Here  the  amplitude  A   (9)   is  the  maximum  value  of  U  ,   the  build-up  time  T     (e)   is 
the  time  at  which  the  maximum  occurs  and  S(T/T      )   is  the  pulse  shape   factor,     A, 


(28)       A. 


(29)       T^„  =  —^ e^  =    (1.85  X  10"^)  9^^ 


(30) 


S(T/T     )   =   (T     A)3/2  ^^  3[^_(^^)l/2j^ 


A  graph  of  the  pulse  shape  factor  S  is  given  in  Figure  1, 

this  graph  by  multiplying  the  vertical  side  by  A  (6)^  which  is  given  by  (28), 

If  the  observation  point  is  above  the  ground  (r  >  a),  U     is  obtained  by 
multiplying  (2?)  by  the  time  dependent  height-gain  factor  H  (tA  -|(r))  defined  by 

v[T^^(r)A-2^/3|ycO)j 
el^^^  "~JI?73pp 


^''^       "^^  ^^  '  .Lol/3,roo 


In  (31)  the  time  T^^^Cr)  is  defined  by 

(32)  T^^ir)  -  2^/3|r^|Q(r-a)/3c. 

A  graph  of  H  (t/T  ,(r))  is  shown  in  Figure  2.     If  the  source  is  above  the  grovmd 
(p  >  a)  a  similar  height-gain  factor  must  be  introduced.     Thus  in  general  we  have 
for  the  perfectly  conducting  case 

(33)  U^(r,e,t)  =  A^(©)S(T/T^^)Hg(TA^^t))Hg(TAg^(p)). 

A  graph  of  U     is  shown  in  Figure  3. 

If  the  times  T  ,(r)  and  T  -(p)  are  both  small  compared  to  the  build  up 

time  T     ,  both  height-gain  factors  are  effectively  equal  to  unity.     Then  U     is 
eo  e 

essentially  the  same  as  it  is  for  both  source  and  receiver  on  the  ground.     On 
the  other  hand  if  T  ,(r)  and  T  ,(p)  are  large  compared  to  T     ,  the  maximum  of 
U     is  less  than  its  value  for  r  -  a,   p  =  a  by  the  factor  H^(T^A„T(r))H  (T  A^(p), 
In  this  case  the  shape  of  the  pulse  is  also  slightly  altered. 


In  case  of  finite  conductivity,  (23)  shows  that  U  is  obtained  by  multi- 
plying (33)  by  the  conductivity  factor  Z^ii:/!:^^,^/!:^) .     This  factor  is  given  by 


(3U)       C^  -  exp 


3(T^/T) 


1/2 


(j"t*2lr-Vla:iJ 


1  + 


2\t^\'\tl\      2|d|6:|2 


The  quantity  |6  1,  which  depends  upon  T/T^  and  upon  en.,  is  given  by  (2h)  with 
der 

(35)     U  (r,9,t)  =  Uf^(r,9,t)  C^(t/t  ,T/To). 


5=0.  Thus  we  have,  denoting  by  U   the  result  (33)  for  perfect  conductivity, 


e'^'^eo'^'"2' 


If  we  insert  the  expi^ssion  (33)  for  U   this  becomes 


(36)  Ug(r,&,t)  -  Ag(©)s(T/Tg^)  Hg(TAg^(r))  Hg(T/Tg^(p))  Cg(T/Tgo,  T/Tg). 

A  graph  of  C  is  shown  in  Figure  U  and  one  of  U  in  Figure  5. 

^  quite  similar  analysis  is  also  possible  in  the  magnetic  case.  From  (22) 
we  find  that  the  analogue  of  (33)  for  U  in  the  case  of  perfect  conductivity  is 


(37)   U^(r,9,t)  »  A^(9)  S(T/T^^)  H^(TA^(r))  HJt/T^(p)). 


The  shape  factor  3  is  the  same  as  in  the  electric  case  but  the  amplitude  A   (  0 ) , 

the  build-up  time  T  ^,   the  times  T  ,(r)   and  T  ,(p)   and  the  height-gain  factors  H 
mo  mj.  mj.  in 

are  slightly  different.     These  quantities  are  given  by 


(38)     A_(e) 


3^0 


"'  '    oii/2  oi/5  ,ji^o,L  _2  •  ,rr 


2^'^  2^'^  ne^lt 


fi 


sin  9 


(3.78x3D"^)       ^        -  (3.78xlC>-3)^, 


^I'^ol^^    3  Q       . 

^^^^       ^mo  "     r ®     "    ^^-^^  ^  '^^     ">  ^; 

3  c 
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(UO)       T  ,(r)   -  T^'^XZlX  9(r-a)/3c, 


v[T^(r)/T-2^/3|r°|] 


(Ul)        H(T/T.(r))=      L"^..     ^., --- 


A  graph  of  the  height-gain  factor  H^  is  given  in  Figure  6.     This  height-gain 
factor  is  zero  at  the  ground.     Therefore  graphs     of  U     ai^  given  for  r  >  a  and 
p  >  a  in  Figure  7. 

In  the  case  of  finite  conductivity,  we  find  from  (2$)  that  U     can  be 
obtained  by  multiplying  the  result   (37)  by  the  conductivity  factor  C   (t/T     .T/Tp). 


This  factor  is 

3xp 
C  li?  ■)       0  C  T  /t      t  A'  ^= 

2i':i'(it° 

-16^1)  ♦  1 

The  quantity   |6   |   depends  upon  tAo  and  upon  O",  .     It  is  given  by  (26)  with  s  »  0. 
If  we  denote  the  result  (37)  for  perfect  conductivity  by  U       we  now  have 

(U3)       Uj^(r,9,t)   =  U^(r,e,t)  CJt/T^,  iH^. 

If  we  replace  U°°  by  means  of  (37)  this  becomes 

(Ul.)       Ujr,e,t)  .=  A^(e)  3(T/T^^)   H^(T/V(r))H^(T/T„l(p))cjT/T^^,T/T2). 


If  o~=  U  mhos/meter,  which  is  the  conductivity  of  sea  water,  C  differs  from 
unity  by  less  than  one  per  cent  for  all  values  of  9.  Therefore  in  this  case  U 
is  essentially  the  same  as  0  ,  which  is  shown  in  Figure  1. 

The  field  components  can  be  found  from  the  Hertz  vector  by  equation  (l)  in  the 
electric  case  and  by  equation  (21)  in  the  magnetic  case.  In  the  electric  case  E^sO, 
and  Eq  "  0  at  r  '  a  in  the  perfectly  conducting  case.  The  main  non-zero  component 
of  E  is  E  .  By  using  (l)  and  (18)  we  find,  for  finite  or  infinite  conductivity,  that 


(U5)  E,.-i 


^1.. 


11 


In  the  maf;netic  case  E     =  E_  e  0, 
only  non-zero  component  of  E, 


From  (2)   and  (l8)  we  find   for  the 


(Ii6) 


=  ii£ 
a 


K\^2 


When  the  source  has  a  time  dependence  f(t)  instead  of  6(t),  the  result  for 
the  Hertz  vector  or  the  field  components  can  be  obtained  from  those  for  the  delta 
function  source.  If  U  denotes  the  Hertz  vector  (electric  or  magnetic)  due  to 
the  source  f,  and  U  the  corresponding  Hertz  vector  due  to  the  delta-function 
source,  then  by  superposition 

00 


(U7) 


U^(r,9,t)  "    f(r)  U^(r,9,t-r)dr. 


From  this  equation  we  see  that  the  minijmam  resolution  time  in  U  is  the 

appropriate  build-up  time  T   or  T  .  Thus  frequencies  higher  than  lA   or 
eo     mo  eo 

lA   are  essentially  lost  from  the  pulse.  We  may  describe  this  by  saying  that 
diffraction  has  the  effect  of  a  low  pass  filter.  This  is  understandable  since  high 
frequency  fields  do  not  so  readily  diffract  around  the  earth. 

An  alternative  to  (U7)j  which  may  be  easier  to  compute  from,  is 


(U8) 


l/(r,e,t)  -  7(ck*)  U^  (r,9,t). 


Here  7  is  the  Fourier  transform  of  f(t)  and  k  is  the  value  of  k  at  the  saddle 
point.  It  is  given  by 

3/2 


(U9) 


ck 


ini  T, 


It 


Here  Tl  is  2^°°  in  the  electric  case  and  t?  In   the  magnetic  case.  Although  (U?) 
0      o  0 

is  valid  for  arbitrary  f(t),  (U8)  applies  only  when  f(t)  represents  a  pulse.  This 
limitation  results  because  the  phase  of  T  was  ignored  in  determining  the  saddle  point, 

Finally  we  must  point  out  a  limitation  on  our  result  (33)  in  the  electric 
dipole  case  for  a  finitely  conducting  earth.  This  limitation  pertains  to  small 
values  of  T.  For  such  values  of  T  the  stationary  value  of  k  or  <a  is  large  in  the 
integral  (18)  for  U  .  But  our  approximate  calculation  of  f  ,  which  occurs  in 
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that  integral,   is  based  on  an  expansion  for  large  cr   and  it  is  not  valid  when  co 
is  also  large.     Therefore  the  result  (23)   is  not  accurate  for  small  T.     However 
this  range  of  T  is  confined  to  a  small  interval  which  terminates  long  before 
the  maximtun  occurs,   in  the  case  we  have  considered.     No  such  limitation  occurs 
in  the  magnetic  dipole  case,  however,   since  then  the  expansion  for  large  <r 
becomes  more  accurate  when  »  is  large. 
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